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Particle Trajectory and Icing Analysis of the E 3 Turbofan Engine 

Using LEWICE3D Version 3 

Colin S. Bidwell 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

Particle trajectory and ice shape calculations were made for the Energy Efficient Engine (E 3 ) using 
the LEWICE3D Version 3 software. The particle trajectory and icing computations were performed using 
the new “block-to-block” collection efficiency method which has been incorporated into the LEWICE3D 
Version 3 software. The E 3 was developed by NASA and GE in the early 1980’s as a technology 
demonstrator and is representative of a modern high bypass turbofan engine. The E 3 flow field was 
calculated using the NASA Glenn ADPAC turbomachinery flow solver. Computations were performed 
for the low pressure compressor of the E 3 for a Mach 0.8 cruise condition at 1 1,887 m assuming a 
standard warm day for three drop sizes and two drop distributions typically used in aircraft design and 
certification. Particle trajectory computations were made for water drop sizes of 5, 20, and 100 pm. 
Particle trajectory and ice shape predictions were made for a 20 pm Langmuir-D distribution and for a 
92 pm Super-cooled Large Droplet (SLD) distribution with and without splashing effects for a Liquid 
Water Content (LWC) of 0.3 g/m 3 and an icing time of 30 min. The E 3 fan and spinner combination 
proved to be an effective ice removal mechanism as they removed greater than 36 percent of the mass 
entering the inlet for the icing cases. The maximum free stream catch fraction for the fan and spinner 
combination was 0.60 while that on the elements downstream of the fan was 0.03. The non-splashing 
trajectory and collection efficiency results showed that as drop size increased impingement rates 
increased on the spinner and fan leaving less mass to impinge on downstream components. The SLD 
splashing case yielded more mass downstream of the fan than the SLD non-splashing case due to mass 
being splashed from the upstream inlet lip, spinner and fan components. The ice shapes generated 
downstream of the fan were either small or nonexistent due to the small available mass and evaporation 
except for the 92 pm SLD splashing case. Relatively large ice shapes were predicted for internal guide 
vane #1 and rotor #1 for the 92 pm SLD splashing case due to re-impingement of splashed mass. 

Introduction 

Over the last several years work has been ongoing to develop tools to analyze aircraft configurations 
subjected to Super-cooled Large Droplet (SLD) conditions (Ref. 1) and High Ice Water Content (HIWC) 
conditions (Refs. 2 to 3). Both the SLD and HIWC environments contain conditions outside of the FAA 
Appendix C (Ref. 4) Certification Envelope. The SLD environment contains larger drops than the 
Appendix C envelope which impinge further back and which are subject to breakup and splash at impact. 
The HIWC environment contains large ice or mixed phase particles (>50 pm) in very high concentrations 
(~10 g/m 3 ) up to very high altitudes (-40000 ft). This HIWC environment has been responsible for many 
engine anomalies including engine roll backs and shutdowns. Work is underway to quantify the HIWC 
environment and to develop ground test facilities and computational tools to assess the sensitivity of 
various engines to the HIWC threat. New certification rules, which will require aircraft to fly safely 
through these conditions, are on the horizon. These new rules will generate new requirements for aircraft 
ice protection system design and certification and for the tools which aid in these processes. 
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The development of icing computational analysis tools which produce sufficiently accurate results in 
a reasonable amount of computational time for turbomachinery flows has been a major challenge 
(Refs. 5 and 6). The use of unsteady tools to simulate the flow and particle transport in the highly time 
dependent turbomachinery flows was seen as impractical and as possibly unnecessary. For this reason a 
methodology was developed at NASA Glenn which uses the steady flow assumption commonly used in 
turbomachinery design tools. These methods typically model blade rows as a single blade with 
circumferential symmetry and circumferentially averaged inflow and outflow boundary conditions which 
are generated from neighboring blade rows. These methods typically march through the turbomachinery 
calculating steady flow for each blade using the upstream blade outflow boundary data and the 
downstream inflow boundary data for the inflow and outflow boundary conditions respectively. Several 
passes through the engine to achieve convergence are typically employed by these methods. 

The newly developed NASA Glenn “block-to-block” icing analysis method follows the same logic 
used in these steady flow turbomachinery flow analysis tools. Droplet transport and ice shape calculations 
are generated for each blade row using steady, single blade flow solutions and the outflow droplet size, 
concentration, and particle velocity data from the upstream blade row as inflow data. The upstream blade 
outflow droplet concentration and droplet velocities are circumferentially averaged before being passed to 
the downstream blade row as an inflow boundary condition. This “block-to-block” method was 
incorporated into the NASA Glenn 3D ice accretion tool LEWICE3D Version 3 (Ref. 7). 

The E 3 engine (Refs. 8 to 10) was selected as a test case for the newly developed “block-to-block” 
method incorporated into LEWICE3D Version 3. The E 3 was developed by NASA and GE in the early 
1980’s as a technology demonstrator. The engine was chosen because it is representative of a modem 
high bypass turbofan engine, the geometry and experimental data were publicly available and flow 
solutions were readily available. 

The analysis was done in two parts. The first part was a particle analysis for 5, 20, and 100 pm 
particles to explore the nature of particle flow through the E 3 . The particles used in this analysis were 
water droplets but they also could represent spherical ice particles which have a similar density. 
Throughout the study the term particle and drop are used interchangeably and are meant to imply particles 
with the same drag characteristics as water drops with a density equal to that of water. The 5 and 20 pm 
particles are typical of particle sizes used in Appendix C analysis. The 100 pm particle is a typical particle 
size for both SLD and HIWC. The second part was an icing analysis using two water droplet distributions 
typically used in aircraft design and certification. The first distribution was a 20 pm Appendix C 
distribution. The second distribution used was a 92 pm SLD distribution. The SLD droplet distribution 
was obtained from measurements in the NASA Glenn Icing Research Tunnel during droplet impingement 
tests (Ref. 11). The SLD calculations were done with the SLD splashing model turned off and with it 
turned on to quantify the effect of droplet splashing on the icing characteristics of the E 3 . The SLD 
splashing results may provide some insight into the HIWC ice crystal environment throughout the engine. 
It is thought that the large ice crystals contained in the HIWC environment probably breakup into smaller 
particles during impact much like the large SLD drops breakup into smaller drops upon impact. It is also 
important to note that the E 3 flight condition chosen for the analysis (cruise at Mach 0.8, 1 1,887 m) is 
typical of HIWC conditions but not of SLD or Appendix C icing conditions, which typically occur at 
lower altitudes and speeds. The authors believe the icing analysis has merit even though the conditions 
chosen were not consistent with actual Appendix C or SLD conditions. The study represents an initial 
step in exploring the particle and icing characteristics of a turbofan compressor and it illustrates the newly 
developed capability of the LEWICE3D Version 3 software. 

Nomenclature 


d droplet diameter, pm 

E 3 Energy Efficient Engine 

HTC Heat Transfer Coefficient, W/m 3 /K 
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IGV Inlet Guide Vane 
LWC Liquid Water Content, g/m 3 
MVD Median Volume Diameter, pm 
P Collection efficiency 


Analytical Method 

Grid and Flow Calculations 

The ADPAC flow solver (Refs. 12 to 14) was used to generate the flow solution for the analysis. The 
ADPAC code is a three-dimensional, finite volume based, Reynolds Averaged Navier-Stokes flow solver. 
The code computes flows on complex propulsion system configurations using multi-block body fitted 
grids. The method employs a “mixing-plane” procedure to pass boundary condition data between grid 
blocks for the steady state flows. The code supports parallel computing and uses a Baldwin-Lomax based 
turbulence model. 


LEWICE3D 

The LEWICE3D ice accretion code was used for the drop trajectory and icing analysis. The grid 
based code incorporates droplet trajectory, heat transfer and ice shape calculation into a single computer 
program. The code can handle generic multi-block structured grid based flow solutions, unstructured grid 
based flow solutions, simple Cartesian grids with surface patches, and adaptive grids with surface 
patches. The latter two methods allow the use of generic panel code input which, when combined with 
LEWICE3D, is a computationally efficient method for generating ice shapes. The code can handle 
overlapping and internal grids and can handle multiple planes of symmetry. Calculations of arbitrary 
streamlines and trajectories are possible. The code has the capability to calculate tangent trajectories and 
impingement efficiencies for single drops or drop distributions using area based collection efficiency 
methods. Ice accretions can be calculated at arbitrary regions of interest in either a surface normal or 
tangent droplet trajectory direction. The program can run on a variety of single processor and parallel 
computers, including Unix, Linux, and Windows (Microsoft Corporation) based systems. 

Version 3 of the LEWICE3D software, which incorporates several new features, was used for the 
analysis. These features include a new particle splash and bounce algorithm, a new geometry handling 
scheme which allows complex mirroring, transformation and relative motion of input grid blocks and a 
new algorithm which calculates block-to-block collection efficiencies. The capability to calculate particle 
trajectories in rotating reference frames was also added (Ref. 5). These new additions will enable users to 
analyze SLD conditions and to calculate collection efficiency with particle splash and bounce effects 
through turbomachinery. 

The “block-to-block” collection efficiency methodology in LEWICE3D was used to calculate 
collection efficiencies and ice shapes for the low pressure compressor of the E 3 for a cruise condition 
(Mach 0.8, at an altitude of 1 1,887 m assuming a standard warm day). The flow solution was obtained 
using the NASA Glenn ADPAC flow solver. The work presented in this paper represents calculations on 
the E 3 low pressure compressor for several drop sizes and droplet distributions. 

Analysis 

The E 3 analysis included the calculation of flow, particle trajectories and ice accretions. The results 
for the flow are presented along with particle analysis for 5, 20, and 100 pm drop sizes. Ice accretion 
calculations were made for 2 droplet distributions typically used in aircraft design and certification. These 
distributions included an Appendix C based 7 bin Langmuir-D distribution with a Median Volume 
Diameter (MVD) of 20 pm (Table 1) and a 10 bin SLD based distribution with an MVD of 92 pm 
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(Table 2). Ice accretion calculations were made for the SLD distribution with and without splashing. The 
ice accretion calculations were made at six locations in the E 3 low pressure compressor. These included 
the inlet lip, splitter lip #1, splitter lip #2, and the 25 percent span location on the fan blade, inlet guide 
vane #1, and rotor #1 shown in Figure 1. 

The grid structure used for the flow and particle analysis is shown in Figure 1. The grid contained 
12 structured, abutted grid blocks with a total of 327,583 nodes. The steady, viscous flow solution was 
generated for a Mach 0.8 cruise condition at 1 1,887 m assuming a standard warm day. Flow vectors along 
the centerline of the axi-symmetric solution are shown in Figure 2. 

Figure 3 shows the collection efficiency (|3) and particle trajectory results for the 5 pm case. 

C ollection efficiency is a non-dimensional measure of the water flux for a surface and is dependent upon 
the amount of convergence or dispersion of particles in a flow and the orientation of the surface relative to 
the droplet paths. A collection efficiency of one implies the surface water loading is equal to the free 
stream water loading. A value less than one means the surface water loading is less than the free stream 
water loading and a value greater than one imply that the surface water loading is greater than the free 
stream level. From the collection efficiency, drop trajectory and drop trajectory impact points shown in 
Figures 3(a) to (d) we can see that impingement occurs throughout the low pressure compressor. The 
droplet trajectories are shown in the rotating reference frame of the fan. The droplet trajectory impact 
point plot (3d) displays the impact locations for the droplets calculated in each of the blocks for the 
“block-to-block” method. This plot illustrates the density and location of the impacting drops used in the 
“block-to-block” method. The peak value of collection efficiency for the fan was 1.2 (Table 3) and 
occurred on leading edge of the pressure side of the fan blade near the blade root. The free stream catch 
fractions were relatively small for the 5 pm for the various components shown in Table 3 (<6 percent). 
The free stream catch fraction was defined as the ratio of the mass impinging on a component divided by 
the mass available in the free stream for an area equal to the area bounded by the highlight of the inlet lip. 
The elements of interest listed in Table 3 (spinner, fan blade, splitter lip #1, IGV #1, rotor #1 and splitter 
lip #2) removed 1 1 percent of the mass entering the inlet with the spinner and fan removing the largest 
portion (7 percent). 

The collection and particle trajectory results for the 20 pm case are show in Figure 4. As for the 5 pm 
drop the 20 pm drop revealed impingement throughout the low pressure compressor although the 
impingement area was larger for the 20 pm drop (Figs. 3(d) and 4(d)). For most of the components the 
impingement rates were higher for the 20 pm drop than for the 5 pm drop due to drop inertia. Larger 
particles, which have larger inertia, are more resistive to changes in direction due to flow gradients than 
smaller particles which results in the larger particles being less apt to avoid obstacles. The maximum 
value of collection efficiency for the fan was 1.54 (Table 4) and it occurred on leading edge of the 
pressure side of the fan blade near the root of the blade. The elements of interest listed in Table 4 
(spinner, fan blade, splitter lip #1, IGV #1, rotor #1 and splitter lip #2) removed 49 percent of the mass 
entering the inlet. The fan and spinner elements removed most of this mass (47 percent). 

The collection and particle trajectory results for the 100 pm case are show in Figure 5. Impingement 
is visible from Figure 5(d) for the bypass region but not aft of IGV #1 for the inner flow path. The 
100 pm drops were not able to negotiate the flow path between the IGV #1 blades due to droplet inertia 
and a higher relative angle-of-attack. From the IGV #1 inflow particle and air velocity vectors shown in 
Figure 6 we can see that that the 100 pm drop had a larger relative angle-of-attack than 5 pm particle or 
the air. This higher relative angle-of-attack and inertia of the larger 100 pm drops results in a 100 percent 
collision rate with IGV #1. The maximum collection efficiency on the fan blade was 0.92 and occurred on 
the aft end of the fan blade near the root. The elements of interest listed in Table 5 (spinner, fan blade, 
splitter lip #1, IGV #1, rotor #1 and splitter lip #2) removed 73 percent of the mass entering the inlet with 
the spinner and fan removing the largest portion (72 percent). 

The collection efficiency results for the 20 pm Langmuir-D distribution are shown in Figure 7. The 
results are similar to those for the single 20 pm particle. The impingement limits are farther aft on the 
inlet lip and spinner for the distribution than for the single drop size case due to the presence of larger 
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drops in the distribution. The maximum collection efficiency on the fan blade (Table 6) was also greater 
for the distribution case than that for the single drop size case due to the presence of larger drops in the 
distribution case (1.9, 1.5 respectively). The maximum occurred on the leading edge of the pressure side 
of the fan blade near the blade root. The amount of mass impinging on the elements of interest in Table 6 
was also similar to the 20 pm single drop case (47 and 49 percent, respectively). The distribution case 
yielded a small amount of mass entering the inner core (0.2 percent of that entering the inlet). This was 
due to the presence of small drops in the distribution case which were able to negotiate the blade passages 
close to the inner wall and enter the inner core. 

Figure 8 presents the collection efficiency results for the 92 pm SLD distribution without splashing. 
The magnitude of collection efficiency and the extent of impingement on the inlet lip and spinner are 
larger than for the 20 pm Langmuir-D case due to the larger drops present in the SLD distribution. The 
impingement rates and maximum collection efficiencies for the 92 pm distribution case were similar to 
those for the single 1 00 pm drop case except that the distribution case generated impingement on the 
elements aft of IGV #1. This was due to the presence of smaller drops in the distribution case which were 
able to negotiate the passages between the upstream blades. The maximum collection efficiency on the 
fan blade was 1.14 and occurred on the aft end of the fan blade near the root. The elements of interest 
listed in Table 7 (spinner, fan blade, splitter lip #1, IGV #1, rotor #1 and splitter lip #2) removed 
67 percent of the mass entering the inlet with the spinner and fan removing the largest portion 
(66 percent). As for the 20 pm Langmuir-D distribution case a small amount of mass was calculated 
entering the inner core (0.03 percent of that entering the inlet). 

The SLD 92 pm collection efficiency results for splashing are shown in Figure 9 and in Table 8. The 
SLD splashing model results in an increase in mass loss as the impingement limits are approached for the 
inlet lip and spinner over the non-splashing result shown in Figure 8. The fan blade impingement 
characteristics for the SLD splashing and non splashing cases were markedly different. The splashing 
case shows a large amount of mass loss near the leading edge of the fan blade toward the root. The 
splashing case also shows higher levels of collection efficiency on the trailing edge of the fan blade near 
the root than the non-splashing case due to re-impingement. The mass collected on the fan blade is 
smaller for the SLD splashing case than for the non-splashing case due to splashing (91.466 and 
134.717 g/s, respectively). The percentage of mass entering the inlet which was collected on the elements 
of interest listed in Table 8 (spinner, fan blade, splitter lip #1, IGV #1, rotor #1 and splitter lip #2) was 
much smaller for the splashing case than for the non-splashing case due to splashing (39 and 67 percent, 
respectively). The splashing case also predicted more mass entering the core than the non-splashing case 
(6.946 and 0.061 g/s, respectively). This additional mass was comprised of a large number of small 
particles generated by the breakup of the larger SLD drops during impact. 

The ice accretion results for several elements of interest in the E 3 compressor are presented in 
Figures 10 to 23. These results include collection efficiency, heat transfer coefficient and ice shape 
predictions for the inlet lip, fan blade, splitter lip #1, IGV #1, rotor #1, and splitter lip #2. The predictions 
were generated for an LWC of 0.3 g/m 3 , an icing time of 30 min for a 20 pm Langmuir-D distribution, a 
92 pm SLD distribution without splashing and a 92 pm SLD distribution with splashing. 

From the collection efficiency results for the inlet lip in Figure 10 we can see that the SLD droplet 
distributions yield larger maximum collection efficiencies and extents of impingement than the 
Langmuir-D distribution due to the presence of larger drops. The SLD splashing model also generates 
lower impingement efficiencies than the SLD non-splashing model. 

The fan blade collection efficiency distributions are shown in Figure 11. The SLD non-splashing case 
yielded the largest overall collection. The overall collection was defined here as the area under the 
collection efficiency curve. The SLD splashing case resulted in lower collection efficiency at the leading 
edge of the fan blade than for the SLD non-splashing case due to drop splashing and a higher level of 
collection efficiency at the trailing edge due to re-impingement of mass. The SLD cases showed an 
increase in collection efficiency towards the trailing edge of the fan blade while the 20 pm Langmuir-D 
distribution was relatively flat. 
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The collection efficiency distributions for splitter lip #1 are shown in Figure 12. The 20 pm 
Langmuir-D distribution generated the highest collection efficiency on splitter lip #1 although the level 
was small (P <0.1 1). From Tables 6 to 8 we see that the impingement rates were also higher for the 
smaller Langmuir-D drop size case than for the SLD non-splashing and splashing cases (0.645, 0.401, and 
0.105 g/s, respectively). The SLD splashing case revealed a large loss in mass over the non-splashing case 
due to droplet splashing. The smaller impingement rates of the larger 92 pm non-splashing SLD case as 
compared to the smaller 20 pm Langmuir-D case was counter-intuitive. Generally larger drops produce 
larger collection efficiency on static objects. Two contributing factors to this result are postulated. The 
first is that there is less mass available for impingement at the inflow to splitter lip #1 for the larger drops 
than for the smaller drops because the fan and spinner have removed more of the mass for the larger 
drops. The fan and spinner removed 67 percent of the mass entering the inlet for the SLD non-splashing 
case as opposed to 47 percent for the 20 pm Langmuir-D case. It was also thought that the reduction in 
collection efficiency could be due in part to a sweep effect caused by the swirl of the particle paths by the 
fan. From Figure 6 we can see that the larger 100 pm drops have a larger initial cross flow component of 
velocity relative to splitter lip #1 as compared to the smaller 5 pm drops. This larger cross flow angle or 
relative sweep could result in a reduction of collection efficiency as it does for external swept wing 
droplet impingement (Ref. 15). 

The collection efficiency results for IGV #1 are shown in Figure 13 and in Tables 6 to 8. The SLD 
splashing case shows the largest impingement rate due to re-impingement of mass splashed from 
upstream components. Impingement on both the pressure and suction sides of IGV #1 were calculated for 
all of the distribution cases. For the non-splashing cases most of the impingement was on the upper 
suction surface due to the low relative angle-of-attack of the inflow particles (Fig. 6). The 92 pm SLD 
non-splashing case yielded a smaller impingement rate than the 20 pm Langmuir-D case (0.9 and 
1.1 percent of the mass entering the inlet, respectively) because of the smaller available mass due to larger 
impingement on upstream components for the 92 pm SLD case. The large amount of re-impingement on 
the trailing edge of the pressure surface for the SLD splashing case was comprised of the droplets which 
were splashed from the pressure side of the fan blade and hence entered the IGV #1 at a higher relative 
angle-of-attack than for the non-splashed drops. 

Rotor #1 collection efficiency results are shown in Figure 14. For all three cases there was 
impingement on the entire length of the pressure side of rotor #1 and on the leading portion of the suction 
side. The relative amount of impingement was small for the non-splashing cases (|3 <0.06). As for IGV #1 
the 92 pm SLD non-splashing case yielded a smaller impingement rate than the 20 pm Langmuir-D case 
(0.06 and 0.55 percent of the mass entering the inlet, respectively) because of the smaller available mass 
for the 92 pm SLD case. We can see that for the rotor #1 case the SLD splashing case yields 
1165 percent (Tables 7 and 8) more impingement than the SLD non-splashing case due to re- 
impingement. The re-impinging mass for the rotor was a combination of the mass splashed from the 
upstream components including the spinner, the fan, and IGV #1. 

The collection efficiency results for splitter lip #2 are shown in Figure 15 and Tables 6 to 8. The 
maximum collection efficiencies ([1 <0.08) and free stream catch fractions for the splitter #2 lip were 
small for all of the cases (<0.00 13). The 92 pm SLD non-splashing case yielded a smaller impingement 
rate than the 20 pm Langmuir-D case (0.03 and 0.17 percent of the mass entering the inlet, respectively) 
because of the smaller available mass for the 92 pm SLD case. The SLD splashing case yields 
443 percent (Tables 7 and 8) more impingement than the SLD non-splashing case due to re-impingement. 

The heat transfer coefficient distributions for the inlet lip, splitter lip #1 and splitter lip #2 are shown 
in Figure 16. The largest values were observed for splitter lip #1 due to the smaller leading edge radius of 
curvature for this element. The smaller leading edge radius for this element translates into large velocity 
gradients which yield higher heat transfer coefficients than for the other two elements. All three of the 
elements exhibit transition both above and below the stagnation point evidenced by sharp increase in the 
heat transfer coefficient on either side of the stagnation point. The splitter lip #1 heat transfer coefficient 
distribution shows the earliest transition due to its sharper leading edge. 
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Figure 17 shows the results for the fan blade, IGV #1 and rotor #1 . The highest heat transfer 
coefficients are realized by IVG #1 due to its smaller leading edge radius of curvature. Transition occurs 
immediately for all three of the elements due to the relatively small leading edge radius of curvature of 
the elements. 

Ice shape predictions for the inlet lip are shown in Figure 18. The SLD conditions generated the 
largest ice shapes due to the higher collection efficiencies produced by the larger drops present in the 
SLD distributions. The non-splashing SLD case results in a larger shape than the splashing SLD case due 
to mass loss by droplet splashing. The higher water loading for the SLD cases results in a glaze condition 
with small horns evident in the figure. The smaller 20 pm Langmuir-D generates a lower water loading 
and a shape more indicative of a mixed condition. 

The ice shape comparisons for the fan blade are shown in Figure 19. The ice shapes for the non- 
splashing cases were very similar with the smaller 20 pm Langmuir-D case generating a smaller trailing 
edge ice thickness due to lower collection efficiencies than for the SLD non-splashing case. The effect of 
mass loss at the leading edge due to splashing is evident from the splashing and non-splashing ice shapes. 
Also from the figure one can see the increased mass at the trailing edge for the SLD splashing case over 
the SLD non-splashing case due to re-impingement. 

The ice shape predictions for splitter #1 are shown in Figure 20. The largest ice shapes were produced 
by the 20 pm Langmuir-D case which also showed the highest collection efficiencies. The SLD splashing 
case is markedly smaller than the SLD non-splashing case due to drop splashing and the resulting mass loss. 

The ice shapes for IGV #1 are shown in Figure 21. The largest ice shape was produced by the 
splashing SLD shape due to the larger impingement rates for this case. The SLD splashing case generated 
a large amount of trailing edge ice on the pressure side of the blade due to re-impingement. The SLD non- 
splashing produced a small ice accretion due to the low overall collection efficiency for the case and 
evaporation. All of the ice shapes for IGV #1 were affected significantly by evaporation. The SLD non- 
splashing case lost 69 percent while the SLD splashing case lost 1 9 percent of the impinging mass due to 
evaporation. The 20 pm Langmuir-D case lost 48 percent of the mass due to evaporation. 

Ice shape predictions for rotor #1 are shown in Figure 22. The ice shape produced for the SLD 
splashing case was relatively large while those for the SLD non-splashing and 20 pm Langmuir-D cases 
were small. The percent evaporation for the non-splashing SLD case was 100 percent while that for the 
splashing case was 15 percent. The 20 pm Langmuir-D case lost 70 percent of the mass due to 
evaporation. 

The splitter lip #2 ice accretions are shown in Figure 23. Relatively small accretions were generated 
for the SLD splashing and 20 pm Langmuir-D cases. The largest ice shape was produced for the SLD 
splashing case due to re-impingement of mass splashed from upstream components. The SLD non- 
splashing case, which showed some impingement, lost 100 percent of the impinging mass due to 
evaporation. The evaporation for the SLD splashing case was 58 percent while that for the 20 pm 
Langmuir-D case was 79 percent. 


Conclusion 

Predictions for droplet impingement and ice accretion were generated for the E 3 low pressure 
compressor using the new “block-to-block” collection efficiency method incorporated into LEWICE3D 
Version 3 and a flow solution from the ADPAC flow solver. The non-splashing trajectory and collection 
efficiency results showed that as drop size increased impingement rates increased on the spinner and fan 
leaving less mass to impinge on downstream components. The SLD splashing case yielded more mass 
downstream of the fan than the SLD non-splashing case due to mass being splashed from the upstream 
inlet lip, spinner and fan components. The maximum free stream catch fraction for the fan and spinner 
combination was 0.60 while that on the elements downstream of the fan was 0.03. 
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The ice shapes calculated for the E 3 were highly dependent upon particle size, evaporation and 
splashing. The ice accretions for the inlet lip and fan blade were relatively large and were similar in size 
for the Langmuir-D and SLD distributions. The E 3 fan and spinner combination proved to be a good 
collector removing greater than 36 percent of the mass entering the inlet for the distribution cases. The ice 
shapes predicted for the components downstream of the fan were relatively small due to the availability of 
mass and evaporation except for the IGV #1 and rotor #1 SLD splashing cases. The ice shapes or the 
IGV #1 and rotor #1 SLD splashing cases were relatively large due to the re-impingement of mass 
splashed from upstream components. 
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TABLE 1 . — LANGMUIR-D DROPLET 
DISTRIBUTION FOR AN MVP OF 20 |gm 


D, 

% LWC 

jam 


6.2 

0.05 

10.4 

0.1 

14.2 

0.2 

20 

0.3 

27.4 

0.2 

34.8 

0.1 

44.4 

0.05 


TABLE 2.— SLD DROPLET DISTRIBUTION 


FOR AN MVP OF 92 nm 


D, 

% LWC 

(im 


10.96046 

0.05 

27.59405 

0.1 

51.17964 

0.2 

92.18893 

0.3 

136.2497 

0.2 

179.8884 

0.1 

230.625 

0.03 

283.3301 

0.01 

339.5112 

0.005 

391.7745 

0.005 


NASA/TM— 20 12-2 1 7696 


9 





TABLE 3.— E 3 IMPINGEMENT STATISTICS FOR A DROP SIZE OF 5 pm 


Element 

Pmax 

Impingement rate, 
g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 nr) 

l 

237.397 

1 

Inlet lip 

0.300 

2.186 

0.0092 

Inlet capture 

1.248 

189.932 

0.8000 

Spinner 

0.034 

0.338 

0.0014 

Fan blade 

1.200 

13.655 

0.0575 

Splitter lip #1 

0.305 

1.134 

0.0048 

IGV #1 

0.518 

4.757 

0.0200 

Rotor #1 



0.0038 

Splitter lip #2 



0.0000 

Inner core 

0.513 

3.487 

0.0147 


TABLE 4.— E 3 IMPINGEMENT STATISTICS FOR A DROP SIZE OF 20 pm 


Element 

Pmax 

Impingement rate, 
g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 nr) 

l 

237.397 

1 

Inlet lip 

0.671 

25.395 

0.1070 

Inlet capture 

1.057 

190.128 

0.8000 

Spinner 

0.174 

5.110 

0.0210 

Fan blade 

1.540 

83.894 

0.3530 

Splitter lip #1 

0.183 

1.136 

0.0048 

IGV #1 

0.420 

2.286 

0.0096 

Rotor #1 

0.160 

0.612 

0.0026 

Splitter lip #2 

0.006 

0.046 

0.0002 

Inner core 

0 

0 

0 


TABLE 5.— E 3 IMPINGEMENT STATISTICS FOR A DROP SIZE OF 100 nm 


Element 

Pmax 

Impingement rate, 
g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 m 2 ) 

l 

237.397 

1 

Inlet lip 

0.950 

85.570 

0.3600 

Inlet capture 

1.013 

199.623 

0.8410 

Spinner 

0.528 

28.160 

0.1190 

Fan blade 

0.920 

115.134 

0.4850 

Splitter lip #1 

0.070 

0.433 

0.0018 

IGV #1 

0.160 

1.679 

0.0071 

Rotor #1 

0 

0 

0 

Splitter lip #2 

0 

0 

0 

Inner core 

0 

0 

0 
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TABLE 6.— E 3 IMPINGEMENT STATISTICS FOR A 20 nm MVP LANGMUIR-D DISTRIBUTION 


Element 

Pmax 

Impingement rate, 

g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 m 2 ) 

l 

237.397 

1 

Inlet lip 

0.659 

26.574 

0.1119 

Inlet capture 

1.103 

190.289 

0.8016 

Spinner 

0.178 

6.270 

0.0264 

Fan blade 

1.943 

78.901 

0.3324 

Splitter lip #1 

0.164 

0.645 

0.0027 

IGV #1 

0.203 

2.081 

0.0090 

Rotor #1 


1.054 

0.0044 

Splitter lip #2 


0.318 

0.0013 

Inner core 

0.086 

0.432 

0.0018 


TABLE 7.— E 3 IMPINGEMENT STATISTICS FOR A 92 pm MVP SLD DISTRIBUTION WITHOUT SPLASHING 


Element 

Pmax 

Impingement rate, 
g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 m 2 ) 

l 

237.397 

1 

Inlet lip 

0.890 

77.401 

0.3260 

Inlet capture 

1.030 

199.847 

0.8334 

Spinner 

0.405 

24.243 

0.1021 

Fan blade 

1.140 

108.041 

0.4451 

Splitter lip #1 

0.081 

0.401 

0.0017 

IGV #1 

0.142 

1.860 

0.0078 

Rotor #1 

0.027 

0.116 

0.0005 

Splitter lip #2 

0.008 

0.056 

0.0002 

Inner core 

0.016 

0.061 

0.0003 


TABLE 8.— E 3 IMPINGEMENT STATISTICS FOR A 92 pm MVP SLD DISTRIBUTION WITH SPLASHING 


Element 

Pmax 

Impingement rate, 
g/s 

Free stream catch fraction 

Free stream capture tube 
(assuming capture area of 3.26 m 2 ) 

l 

237.397 

1 

Inlet lip 

0.886 

40.904 

0.1723 

Inlet capture 

1.334 

210.935 

0.8885 

Spinner 

1.563 

14.662 

0.0618 

Fan blade 

3.057 

71.090 

0.2995 

Splitter lip #1 

0.055 

0.105 

0.0004 

IGV #1 

1.178 

3.828 

0.0161 

Rotor #1 

2.133 

1.467 

0.0062 

Splitter lip #2 

0.088 

0.304 

0.0013 

Inner core 

1.288 

6.946 

0.0293 
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(a) Surface model. 


(b) Grid block structure. 


z 



Figure 1. — Surface model and grid structure for E 3 flow model. 
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(a) Side view. 



(b) Orthogonal view. 

Figure 2. — Centerline velocity vectors for the E 3 low pressure compressor. 
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(a) Collection efficiency. 
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(b) Drop trajectories (front view). 

Figure 3. — Drop trajectory and collection efficiency results for the E 3 for 5 pm drop. 
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(c) Drop trajectories (axial view). 
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x 




(d) Drop impact locations (axial view). 


Figure 3. — Concluded. 
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(a) Collection efficiency. 
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(b) Drop trajectories (front view). 


Figure 4. — Drop trajectory and collection efficiency results for the E 3 for 20 pm drop. 


NASA/TM— 2012-217696 


16 






1 

i 


BETA 
1.0000 
0.9583 
0.9167 
0.8750 
0.8333 
0.7917 
0.7500 
0.7083 
0.6687 
0.6250 
0.5833 
0.5417 
0.5000 
0.4583 
0.4167 
0.3750 
0.3333 
0.2917 
0.2500 
0.2083 
0.1667 
0.1250 
0.0833 
0 . 04 1 7 
0.0000 


(c) Drop trajectories (axial view). 
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(d) Drop 


impact locations (axial view). 



Figure 4. — Concluded. 
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(b) Collection efficiency. 
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(b) Drop trajectories (front view). 


Figure 5. — Drop trajectory and collection efficiency results for the E 3 for 100 pm drop. 
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(c) Drop trajectories (axial view). 

z 


x 



(d) Drop impact locations (axial view). 


Figure 5. — Concluded. 
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(a) Front view. 



(b) Axial view. 


Figure 7. — Collection efficiency results for the E 3 for a 20 gm MVD Langmuir-D distribution. 
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(a) Front view. 



(b) Side view. 

Figure 8. — Collection efficiency results for the E 3 for a 92 pm MVD SLD distribution with no splashing. 
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(a) Front view. 



(b) Axial view. 


Figure 9. — Collection efficiency results for the E 3 for a 92 urn MVD SLD distribution with splashing. 
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Figure 1 0. — Collection efficiency distribution for the E inlet lip. 
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Figure 1 1 . — Collection efficiency distribution for the E 3 fan blade at 25 percent span. 
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Figure 12. — Collection efficiency distribution for the E 3 splitter lip #1. 
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Figure 13. — Collection efficiency distribution for the E 3 IGV #1 at 25 percent span. 
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Figure 14. — Collection efficiency distribution for the E 3 rotor #1 at 25 percent span. 
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Figure 15. — Collection efficiency distribution for the E 3 splitter lip #2 at 25 percent span. 
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Figure 16. — Heat transfer coefficient distributions for the E 3 . 
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Figure 17. — Heat transfer coefficient distributions for the E 3 blades at 25 percent span. 
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Figure 18. — Ice shapes for the E inlet lip. 



Figure 19. — Ice shapes for the E 3 fan blade. 
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Figure 20. — Ice shapes for the E 3 splitter lip #1 . 
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Figure 21 . — Ice shapes for the E 3 IGV #1 . 
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Figure 22. — Ice shapes for the E 3 rotor #1 . 



Figure 23. — Ice shapes for the E 3 splitter lip #2. 
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efficiency results showed that as drop size increased impingement rates increased on the spinner and fan leaving less mass to impinge on downstream components. The SLD 
splashing case yielded more mass downstream of the fan than the SLD non-splashing case due to mass being splashed from the upstream inlet lip. spinner and fan components. 
The ice shapes generated downstream of the fan were either small or nonexistent due to the small available mass and evaporation except for the 92 pm SLD splashing case. 
Relatively large ice shapes were predicted for internal guide vane #1 and rotor #1 for the 92 pm SLD splashing case due to re-impingement of splashed mass. 
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